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Abstract. The paramagnetic phase of heavy fermion systems is investigated, using a 
non-perturbative local moment approach to the asymmetric periodic Anderson model 
within the framework of dynamical mean field theory. The natural focus is on the 
strong coupling Kondo-lattice regime wherein single-particle spectra, scattering rates, 
d.c. transport and optics are found to exhibit (u> /lol,T /u>l) scaling in terms of a single 
underlying low-energy coherence scale ujl- Dynamics/transport on all relevant (w,T)- 
scales are encompassed, from the low-energy behaviour characteristic of the lattice 
coherent Fermi liquid, through incoherent effective single-impurity physics likewise 
found to arise in the universal scaling regime, to non-universal high-energy scales; and 
which description in turn enables viable quantitative comparison to experiment. 
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1. Introduction. 

Lanthanide based heavy fermion (HF) metals constitute a major, long studied class of 
correlated electron materials [1-8]. Their behaviour is quite distinct from conventional 
clean metals, the basic physics being driven by strong spin-flip scattering from essentially 
localised /-levels, generating the large effective mass and attendant low-energy scale 
indicative of strong interactions. The low-temperature (T) state is a lattice-coherent 
Fermi liquid with well defined quasiparticles and coherently screened /-spins, crossing 
over with increasing T to essentially incoherent screening via independent Kondo 
scattering, before attaining characteristic clean metallic behaviour. Physical properties 
of HF are in consequence typically 'anomalous': e.g. the resistivity p(T) shows a 
strong non-monotonic T-dependence, while optics often exhibit rich structure from 
the microwave to the near infrared, and pronounced evolution on low temperature 
scales [1-8]. 

Theoretical treatments of HF centre on the periodic Anderson model (PAM), in 
which a non-interacting conduction band hybridizes locally with a correlated /-level in 
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each unit cell of the lattice; or on its strong coupling limit, the Kondo lattice model. 
The absence of exact results (save for some in one dimension, see e.g. [9]) has long 
spurred the search for suitable approximation schemes. One such framework, which 
has had a major impact in recent years, is provided by dynamical mean field theory 
(DMFT, for reviews see [10-13]). Formally exact in the large-dimensional limit, the 
self-energy within DMFT becomes momentum independent and hence spatially local, 
but still retains full temporal dynamics; such that all lattice models map onto an effective 
single-impurity model with a self-consistently determined host [10-13]. 

That raises an immediate question, easier asked than answered: to what extent are 
the properties of real HF materials captured within a DMFT approach to the PAM? 
To answer this clearly requires direct quantitative comparsion of theory to experiment. 
And a prerequisite to that in turn is a method to solve the PAM — which DMFT 
does not per se provide. The latter has of course been studied extensively using 
a wide variety of techniques. Full scale numerical methods include the numerical 
renormalization group (NRG) [14, 15], quantum Monte Carlo [16-18] and exact 
diagonalization [19], while theoretical approaches encompass finite-order perturbation 
theory in the interaction U [20,21], iterated perturbation theory [22,23], the lattice non- 
crossing approximation [24,25] and the average t-matrix approximation [26], large- N 
mean-field theory/slave bosons [27-29], the Gutzwiller variational approach [30,31] and 
the recently developed local moment approach [32-34]. All of these methods naturally 
have their own virtues. But most possess significant, well known limitations [2], be 
it the general inability of perturbative approaches (and in practice quantum Monto 
Carlo) to handle strong interactions; failure to recover Fermi liquid behaviour at low- 
energies as arises in NCA-based approaches, restriction to the lowest-energy Fermi liquid 
behaviour as in large- iV/slave boson mean-field theories, finite-size effects limiting exact 
diagonalization, and so on. 

To enable viable comparison to experiment requires an approach that can 
adequately handle all experimentally relevant energy and/or temperature scales in the 
strongly correlated HF regime of primary interest; and indeed ideally also across the full 
spectrum of interaction strengths, such that intermediate valence and related behaviour 
can likewise be treated. One such is employed here, the local moment approach 
(LMA) [32-34]. Via study of the generic asymmetric PAM, our essential aims are (i) to 
provide a many-body description of dynamical and transport properties of paramagnetic 
HF, notably single-particle dynamics, d.c. transport and optical conductivities; as 
considered here, (ii) To make direct quantitative comparison with experiment. That 
is taken up in the following paper where comparison to transport/optical properties of 
CeB 6 , CeAl 3 , YbAl 3 and CeCoIn 5 is made. 

Some remarks on the LMA are apposite at this point since the paper will 
focus mainly on results obtained using the approach, with minimal technical details. 
Intrinsically non-perturbative and as such capable of handling strong interactions, the 
LMA [32-41] introduces the physically intuitive notion of local moments [42] from the 
outset. This leads directly to a 'two-self-energy' description in which, post mean- 
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field level, the key correlated spin-flip dynamics is readily captured; corresponding 
in physical terms to dynamical tunneling between initially degenerate local moment 
configurations, which lifts the erstwhile spin degeneracy and restores the local singlet 
symmetry characteristic of a Fermi liquid state. As with all techniques for lattice 
models within DMFT, the LMA originated in study of the single-impurity Anderson 
model (AIM) [35-41], where results for dynamics are known to give good agreement 
with NRG calculations [37,38,40], and for static magnetic properties with known exact 
results [41]. The approach has recently been developed to encompass the Anderson 
lattice (PAM); initially for the particle-hole symmetric limit [32,33] appropriate to 
the Kondo insulating sub-class of heavy electron materials, where for all interaction 
strengths the system is an 'insulating Fermi liquid' that evolves continuously from 
its simple non-interacting limit of a hybridization-gap insulator [43]. From this a 
rich description of transport and optical properties of Kondo insulators arises [33], 
particularly in strong coupling where the system is characterized by an exponentially 
small indirect gap scale A g , such that dynamics/transport exhibit scaling as functions 
of (u/A g ,T/A g ). Exploiting that scaling enables direct comparison to experiment 
with minimal use of 'bare' material/model parameters; and in particular for three 
classic Kondo insulators Ce 3 Bi 4 Pt 3 , SmB % and YbB 12 , leads to what we regard as 
excellent agreement between theory and experiment on essentially all relevant energy 
and temperature scales [33]. 

The particle-hole symmetric PAM is of course special, confined as it is to the case 
of Kondo insulators. Most recently the LMA has been non-trivially extended to handle 
the generic asymmetric PAM [34] and hence HF metals (with the insulating symmetric 
limit recovered simply as a particular case). Single-particle dynamics at T = were 
considered in [34], with a natural emphasis on the strongly correlated Kondo lattice 
regime of localised /-electrons but general conduction ('c') band filling n c . The problem 
was found to be characterized by a single low-energy coherence scale ool — the precise 
counterpart of the insulating indirect gap scale A g , and likewise exponentially small in 
strong coupling — in terms of which dynamics exhibit one-parameter universal scaling 
as a function of Co = oo/ool, independently of either the interaction strength or local f/c 
hybridization. With increasing Co dynamics cross over from the low-energy quasiparticle 
behaviour required by and symptomatic of the coherent Fermi liquid state, to essentially 
incoherent single-impurity Kondo scaling physics at high-u; — but still in the cu-scaling 
regime and as such incompatible [34] with a two-scale 'exhaustion' scenario [44]. 

In this paper we extend the work of [34] to finite temperature, thereby enabling 
access to d.c. transport and optics. Our primary focus is again the strongly correlated HF 
regime and attendant issues of scaling/universality (that play a key role in comparing 
to experiment), the paper being organised as follows. The model and a bare bones 
description of background theory is introduced in section 2, together with preliminary 
consideration of transport/optics. Results for the thermal evolution of single-particle 
dynamics and scattering rates, and the connection between the two, are given in section 
3. The d.c. resistivity is considered in section 4, with particular emphasis in this 
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context on the crossover from the low-T = T/ujl coherent Fermi liquid to the high-T 
incoherent regime, and explicit connection to single-impurity scaling behaviour. Optical 
conductivities on all relevant u- and T-scales are investigated in section 5; and the paper 
concludes with a brief summary. 

2. Model and theory 

The Hamiltonian for the PAM is given by H = H c + Hf + H hyb : 

H = tcY< c tc lCT - 1 c lc la + J2{ e f + Vi-Ji-^) flu 

+ Vj2(flc ia + h.c.) (2.1) 

The first two terms represent the uncorrelated conduction (c) band, 
H c (= Xlk o-( e c+ e k) c ko- c ko-)5 w hh c-orbital site energies e c and nearest-neighbour hopping 
matrix element Uj = t, rescaled as t oc t*/\/~Z~ c in the large dimensional limit where the 
coordination number Z c — > oo [10-13]. The third term describes the correlated /-levels, 
Hf, with site energies e/ and on-site Coulomb repulsion U; while the final term Hh y b 
hybridizes the c- and /-levels locally via the matrix element V, rendering the otherwise 
localised /-electrons itinerant. 

The model is thus characterized by four independent dimensionless parameters, 
€c/t*,V/t*,U/t* and e//i* (t* sets the scale for the width of the free conduction band 
and is taken as the basic unit of energy, £* = 1). An equivalent and somewhat more 
convenient set of 'bare' /material parameters (with = 1) is e c ,V,U and rj, where 
r] — 1 + 2ef/U. This parameter space is large, and as such encompasses a wide range 
of physical behaviour for the paramagnetic phases we consider. The system is of course 
generically metallic, with non-integral /-level and c-band occupancies (rif — Y2a(fLfia) 
and n c = ^2 a (Ci a c ia ) respectively). That in turn extends from the trivial case of 
weakly correlated, perturbative behaviour, through intermediate valence to the strongly 
correlated heavy fermion (HF) regime. It is naturally the latter, characterized by a low- 
energy coherence scale u L , that is of primary interest. The HF (or Kondo lattice) regime 
corresponds to essentially localised /-electrons, rif — > 1, but with arbitrary conduction 
band filling n c , the latter being controlled by e c (which determines the centre of gravity of 
the free (V — 0) conduction band relative to the Fermi level). It arises when tf = — \ej\, 
for |e/|/A > 1 and (U - |e/|)/A > 1 (whence -1 < r] < 1); where A = nV 2 d c (0), 
with g?q(cj) the free conduction electron density of states as specified below and u = 
the Fermi level. The heavy fermion regime forms our main focus here; intermediate 
valence behaviour will be discussed in an experimental context in the following paper. 

The exception to the above behaviour arises when rif + n c = 2. Here the system 
is generically a Kondo insulator (see eg [34]), with an indirect gap in both its T = 
single-particle spectrum and optical conductivity [33]; the canonical example being the 
particle-hole symmetric PAM with e c = and e/ = —U/2, where rif — 1 = n c for all 
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U. Just like its metallic counterpart arising for rif + n c ^ 2, the Kondo insulator is 
however a Fermi liquid, evolving continuously with increasing interaction strength from 
its non- interacting limit (in this case a 'hybridization gap insulator' [43]). As such, the 
Kondo insulating state is obtained simply as a particular limit of the underlying theory. 



2.1. Background theory 

A knowledge of local single-particle dynamics and their thermal evolution is well 
known to be sufficient within DMFT [10-13] to determine transport properties (see 
section 2.2 below). Our initial focus is thus on the local retarded Green functions 
G^iu) (<r-> — W(t)({f ia (t), f} ff })) and likewise G^(u) for the c-levels, with corresponding 
spectra D^(u) = --ImG^w) (and v — c or /). 

Some brief comments on the free conduction band are first required (V = 
in equation (|2.1|) . where the c- and /-subsystems decouple); specified by the local 
propagator Qq{uj) with corresponding density of states (dos) d^u). This is given by 

g c (uj) = H(u;+ - e c ) (2.2a) 

= -+ 1 -^T~\ ^ 

lo+ - e c - b {cj) 

with uj + = uj + i0 + , where for arbitrary complex z 

H{z)=f + °°de^ (2.3) 

denotes the Hilbert transform with respect to Po( e )i such that from equation (|2.2qj) . 
d^(oo) = Pq{oj — e c ) corresponds simply to a rigid shift of Pq{uj) by e c . Equation (|2.26j) 
defines the Feenberg self-energy S (u>) [45,46] as used below, with S (ou) = 5[^q] alone 
(since g$ = H(S + l/g^) from equations (2.2)). The free conduction band is thus 
determined by the non- interacting dos po(e) which, modulo the rigid e c -shift, reflects 
the underlying host bandstructure, Po( e ) = ^ 1 J2k^( e ~ e k)- While the formalism 
below holds for an arbitrary po(e), explicit results will later be given for the hypercubic 
lattice (HCL), for which within DMFT [10-13] p (e) = vr" 1 / 2 exp (-e 2 ) is an unbounded 
Gaussian; and the Bethe lattice (BL), with compact spectrum po(e) = (2/tt) (1 — e 2 ) 1 ^ 2 
for |e| < 1 [10-13]. The HCL will in fact be the primary case, because the Bloch states 
characteristic of it ultimately underlie the lattice coherence inherent to low-temperature 
metallic HF behaviour. 

The major simplifying feature of DMFT is that the self-energy becomes momentum- 
independent and hence site-diagonal [10-13]; and since we are interested in the 
homogeneous paramagnetic phase, the local Green functions G^(uj) = G u (u) [y = 
c, /) are also site- independent. Straightforward application of Feenberg renormalized 
perturbation theory [45,46], then gives the G v (uS) as 

G» = ~ ( ^ v2 ~ (2.4a) 

^ + - e/ - E f {u; T) - ^ + _ ee _ s(h>) 
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v , / , ! + -; V T 1 t^ G%u) (2.4c) 

uj^ — tf — 2-if{uj] T) [ u> + — tf — hf{u>;T) 

where £/(u;; T) is the retarded /-electron self-energy (£/(cu; T) = (uj; T) — iXy(co>; T) 
such that £j(c<j;T) > 0). In equations (2.4), S(uj) is the Feenberg self-energy for the 
fully interacting case, with S(uj) = S[G C ] the same functional of G c (uj) as it is of g$ in 
the V = limit. In consequence, G c (uj) is given using equations (j2.4a|) . (2.2), ()2.Hj) as 

G c (uj) = H(i) (2.5a) 

where j(w] T) (= 7^(0;; T) + z7/(o;; T)) is given by 

7 (w; T) = W + - e c - — ^ . (2.56) 

w+-e / -S / (a;;T) 

Let us first point up the physical interpretation of equations (2.4), (2.5). G c (uj) is 
a local propagator, and as such familiarly expressed as G c (uj) = N^ 1 X!k^ c ( e k> with 

the ek-resolved conduction electron propagator G c (e-^; uj) = [uj + — e c — — £ c (c<j; T)] _1 

and the usual conduction electron self-energy £ c (cj;T) thus defined. Since po(e) = 
N^ 1 <K e ~~ e k); it follows directly that 

/+00 
dep (e)G c (e;uj) = (G c (e;uj)) e (2.6) 
-00 

with (G c (e\ L = e; uj) =) G c (e;uj) = [uj + — e c — £ c (co>;T) — e] _1 . But equation ()2.6|) is 
precisely the form equation (|2.5ajl (with equation (|2.Hj) for if (7)), showing that 

G c (e;a;) = [7(a;;r)-e]- 1 (2.7) 

with 7(0*; T) related to the conduction electron self-energy by 

1 (uj;T) = uj + -6 c -J: c (uj;T); (2.8) 

and hence (via equation (|2.56jl ) that T, c (uj;T) = V 2 [uj + — ej — S^(u;;T)] _1 in terms of 
the /-electron self-energy alone (because the /-levels alone are correlated). 

For an arbitrary conduction band (specified by po(e)) equations (2.4), (2.5) are 
central; for given the self-energy £/(cl>; T), and hence 7(0;; T) from equation (|2.56|) . G c {uj) 
follows directly from the Hilbert transform equation (|2.5oj) . and G* (uj) in turn from 
equation ()2.4c|) . That statement hides however the truly difficult part of the problem: 
obtaining the self-energy T,f(u;T). This is not merely a calculational issue, e.g. the 
need to solve the problem iteratively and self-consistently (any credible approximation 
to Yif(uj]T) will in general be a functional of self-consistent propagators). It reflects 
by contrast the longstanding problem of obtaining an approximate £/(u;;T) that, 
ideally: (i) Handles non-perturbatively the full range of interaction strengths, from weak 
coupling (itself accessible by perturbation theory or simple variants thereof [20-23]) all 
the way to the strongly correlated Kondo lattice regime that is dominated by spin- 
fluctuation physics and typified by an exponentially small coherence scale ujl. (ii) 
Respects the asymptotic dictates of Fermi liquid behaviour on the lowest energy (uj) 
and/or T scales - on the order of \u\ < ujl itself - yet can also handle the full uj 
and/or T range; including the non-trivial dynamics that arise on energy scales up to 
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many multiples of ul yet which remain universal (and the existence of which we find 
to dominate transport and optics), as well as the non-universal energy scales prescribed 
by the bare material parameters of the problem. 

The success of any theory naturally hinges on the inherent approximation to 
Hf(u;T). In this paper we employ the local moment approach (LMA) [32-34], for 
it is known to satisfy the above desiderata and to our knowledge is currently the 
only theory that does. It is based on an underlying two-self-energy description — a 
natural consequence of the mean-field approach from which it starts, and from which 
the conventional single self-energy £ j follows — together with the concept of symmetry 
restoration that is central to the LMA generally [35-41]. Full details of the LMA for 
the PAM, including discussion of its physical basis and content, are given in [32-34]. 
In particular the generic asymmetric PAM (as considered here) is detailed in [34] for 
T = 0; and extension of it to finite-T, required to consider transport and optics, follows 
the approach of [33] where the particle-hole symmetric PAM appropriate to the case of 
Kondo insulators was considered. For that reason further discussion of the approach 
is omitted here. The reader is instead directed to [33, 34] on the PAM, from which 
appropriate results will be used when required; and to [35-41] for Anderson impurity 
models per se where details of the LMA, including its stengths and limitations in relation 
to other approaches, are fully discussed. 

2.2. Transport and optics 

As mentioned above, a knowledge of single-particle dynamics is sufficient within DMFT 
to determine q = transport properties [10-13]. This arises because the strict absence of 
vertex corrections in the skeleton expansion for the current-current correlation function 
means only the lowest-order conductivity bubble survives [47], and a formal result for 
it is thus readily obtained. Denoting the trace of the conductivity tensor by a(u; T) (| 
of which, denoted by a(uj;T), provides an approximation to the isotropic conductivity 
of a 3-dimensional system), this may be cast in the form 



with cr = 7re ? _ a y ~ merely an overall scale factor (a is the lattice constant and a 
is typically of order pa 10 4 — 10 5 fi -1 cm~ 1 ). The dimensionless dynamical conductivity 
F(u; T) naturally depends on the lattice type, and for a Bloch decomposable lattice 
such as the HCL is given (with t* = 1) by [10-13,33] 



where f(uj) = [e^^ + l] -1 is the Fermi function. Here (as in equation ()2.6j) ). the notation 
{A(e;w)) e = de p (e)A(e; u>) denotes an average with respect to the non-interacting 
conduction band dos po(e); and the spectral density D c (e;uj) = — (l/ir)hn.G c (e; uj) with 
G c (e;u) = [j(uj; T) — e] _1 from equation (|2.7jl . Physically, Ji(uj; T) (= hxvy{uj\T)) 
represents the cu-dependent conduction electron scattering rate (inverse scattering time) 



a(u;T) = o- F(u]T) 



(2.9) 



Fhcl^ T) = - du) X [fM - /(w a + u)] {D c (e; ^)D c (e; + cu)) e 




(2.10) 
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arising from electron interactions, Ji(uj; T) = 1/t(u; T) (= — ImE c (u;; T) from equation 
(|2.8|l ). It is given using equation (|2.56j) by 

1 ^ 2 £}(cu;T) 

in terms of the /-electron self-energy; a knowledge of which thus determines the 
scattering rates (considered explicitly in section 3.1), and in consequence the dynamical 
conductivity equation (|2.1()j) (noting that D c (e; uj) = Ji(uj; T)^ 1 / {\^r{uj] T) — e] 2 + 

[7/(^;T)] 2 )). 

Results for Fhcl(w, T) obtained using the LMA will be considered in sections 4,5. 
Here we simply point out an exact result, not apparently well known, for the weight 
of the Drude peak in the T = conductivity. At T = 0, scattering at the Fermi level 
is absent since the system is a Fermi liquid, i.e. E^(u; = 0;T = 0) =0 and hence 
7/(0; 0) = 0. The leading low-frequency behaviour of (uj\ T) is given by 

Ef(u; 0) ~ Ef(0; 0) - (J - 1)uj (2.12) 

where Z — [1 — (<9£y(u;; 0)/du) u=0 ] _1 is the usual quasiparticle weight /inverse mass 
renormalization; hence (from equation (|2.56j) ) 7_r(0; 0) = — e c + V 2 /e}, where 

e} = e/ + Sf(0;0) (2.13) 

is the renormalized /-level energy. A straightforward evaluation of equation ()2.10j) for 
T = and to — > then shows that Fhcl(w;T = 0) contains a 6(u) Drude 'peak' (as 
it must, reflecting the total absence of Fermi level scattering and a vanishing T = 
resistivity). Denoted by Fr, ru <i e {uj\ 0), it is given explicitly by 

Zef V 2 
F Drude (u;0) = 6(u) z ^ 2 | y2 Po(~e c + — ) (2.14a) 

or equivalently 

F Drude (u;0) = 6(u) — ~7~T" Po(-e c + 4) ( 2 - Ub ) 

e f j 

where e* f = e* f /V 2 and 

u L = ZV 2 . (2.15) 

Equations (2.14) are exact, and bear comment. In the trivial limit V = 
where (equation(2.1)) the /-levels decouple from the conduction band, the total Drude 
weight is naturally d^(uj = 0), the free conduction band dos at the Fermi level (recall 
<1q(u) = Pq{uj — e c )). For any V ^ 0, the Luttinger integral theorem requires 

\[n c + n f )= I C+/f p (e)de + 9(-S}) (2.16) 

J — oo 

(with 9(— e*f) merely the unit step function). This again is an exact result, proven in [34]. 
It holds for any interaction U, reflecting the adiabatic continuity to the non-interacting 
limit that is intrinsic to a Fermi liquid; and shows in general that (any) fixed total 
filling n c + rif determines — e c + 1/e^ entering equations (2.14). Of particular interest is 
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of course the strongly correlated HF regime, where rif — > 1. Here uj^ in equation (|2.15|) 
(= ZV 2 /t* with = 1) is the coherence scale: exponentially small in strong coupling 
(because Z is), it is the single low-energy scale in terms of which all properties of the 
system exhibit universal scaling (as shown in [32-34] and pursued below). In the HF 
regime, n c itself is moreover given (see [34]) by 



\n c = I p (e) de (2.17a) 



showing that e c and n c are in essence synonymous, n c = n c (e c ) being determined by e c 
alone. Conjoining this with equation (j2.16|) gives 

\n f = I eC+1/e ' Po (e) de + 9(-e* f ) (2.176) 

J — e c 

so as nj — > 1 (the HF regime), e*j = ej(e c ) = t}( n c) is also determined by e c alone, 
and is typically of order unity. It is this that determines e°ji entering equation ()2.146|) 
for Fn ru de{oj] 0); showing in turn that the net Drude weight is itself oc ujl = ZV 2 , and 
hence exponentially diminished compared to the free conduction band limit. 

We add that Kondo insulators, arising generically for nf + n c = 2 as mentioned 
earlier, are also encompassed by the above. Using J^depo^e) = 1, the Luttinger 
theorem equation ()2.16|) shows that nf+n c = 2 arises either for = (for an unbounded 
Po{e)) or for — e c + 1/^/ outside the band edges of a compact po(e); such that in either 
case the Drude weight in equations (2.14) vanishes, symptomatic of the vanishing T = 
d.c. conductivity characteristic of the Kondo insulating state. 

Our focus above has naturally been on the canonical case of a Bloch decomposable 
lattice. For a Bethe lattice by contrast, F(u;T) is given [33] by (c/ equation (|2.10Jl ) 

F BL (oo; T) = — f fax [fM - f(ux + oo)} D^D^ + u) (2.18) 

U J-oo 

where D c (uo) (= (D c (e;uj)) e ) is the local conduction band spectrum. In particular the 
d.c. conductivity at T = follows as F BL {0;0) = [D c (0)} 2 ; which, using £}(0;0) = 
together with equations (|2.2aj) . (2.5), (|2.12jl . is given by 

^L(0;0) = [p (-e c + ^)] 2 - (2.19) 

In contrast to equations (2.14) there is thus no Drude 5(a;)-peak and the T = d.c. 
resistivity is in general finite, reflecting of course that the underlying one-particle states 
of the BL are not coherent Bloch states. Hence, aside from the case of Kondo insulators 
where the BL (like the HCL) does capture the vanishing T = d.c. conductivity 
and indirect-gapped optics characteristic of the insulator [33], the 'joint density of 
states' type formula equation ()2.19|) should not be taken seriously when considering 
transport /optics of real materials on sufficiently low T and/or uo scales (as discussed 
further in section 4). 
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3. Single-particle dynamics 

We turn now to LMA results for single-particle dynamics at finite-T. Our natural focus 
will be the strong coupling Kondo lattice regime (where nf — > 1), characterized by 
the low-energy lattice scale u>l = ZV 2 . This scale is of course a complicated function 
of the bare/material parameters, ujl = ujl(€ c , U, V 2 , rj) (detailed LMA results for it 
are given in [34], and NRG results in [15]). That dependence is however a subsidiary 
issue in comparison to the fact that, because ujl becomes exponentially small in strong 
coupling, physical properties exhibit scaling in terms of it; i.e. depend universally on 
uj/ujl, independently of the interaction strength. 

Universality in strong coupling single-particle dynamics at T = has been 
considered in [34] for the generic PAM; the essential findings of which are first reprised 
for use below, (i) Both the c-electron spectrum D c {oj) (= D c {uj)/t 1f with t* = 1) and the 
/-electron spectrum ttAqD^(uj) (with A = 7rV 2 po(— e c ) introduced in section 2), exhibit 
universal scaling as a function of u = uj/ujl in a manner that is independent of both the 
interaction strength U and hybridization matrix element V. (ii) That scaling depends in 
general only on e c (or equivalently the conduction band filling n c , see equation (|2.17aj) ) 
which embodies the conduction band asymmetry; and on rj = 1 —2\ef\/U reflecting the 
/-level asymmetry. More specifically, (iii) in the coherent Fermi liquid regime arising 
for \u\ < 1, the /- scaling spectra depend only on e c and are in fact independent of rj 
as well as U and V. In this low-a) regime the scaling spectra amount in essence to the 
quasiparticle behaviour (equations (3.11) of [34]) required by the asymptotic dictates 
of low-energy Fermi liquid theory, (iv) For \u\ ^> 1 by contrast the /- scaling spectra 
depend on the /-level asymmetry rj (albeit rather weakly), but are now independent of 
e c and indeed also of the lattice type; and the spectrum contains a long, logarithmically 
slowly decaying spectral tail, (v) The latter behaviour, which sets in progressively for 
\u\ > 1, reflects in turn the crossover to incoherent effective single- impurity physics that 
one expects to arise for sufficiently high uj (and/or T): for |<D| 1 the scaling form 
of the /-spectrum is found to be precisely that of an Anderson impurity model (AIM). 
With increasing u, dynamics thus cross over from the low-energy quasiparticle behaviour 
symptomatic of the lattice coherent Fermi liquid state to single-impurity Kondo scaling 
physics at high Q (and that this crossover occurs in a single uj/ujl scaling regime is thus 
incompatible with the occurrence of 'two-scale exhaustion' [44] as explained in [34]). 

Figure Q summarises representative results for T = scaling dynamics (irrelevant 
non-universal energy scales such as U, t* (= 1) or A are of course projected out in 
scaling spectra [32-34]). The main figure shows f/c scaling spectra for the HCL as 
functions of the scaled frequency Co = u/ul, for r/ = with e c = (dashed, and n c — 1) 
and 0.3 (solid, with n c ~ 0.68). The e c = example corresponds to the particle-hole 
(p-h) symmetric Kondo insulator, whose spectra are thus gapped at the Fermi level 
Co = (with ujl = ZV 2 here corresponding to the insulating gap scale [32,33]). For the 
asymmetric conduction band e c = 0.3 by contrast, the gap (which is well developed in 
strong coupling [34]) moves above the Fermi level; and a sharp lattice-Kondo resonance 
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symptomatic of the HF metal, straddling the Fermi level and of width oc ul, takes its 
place in the /-spectra. The inset shows the /-spectra on a much larger Cj scale; displaying 
the e c -independence of the slow logarithmic tails [34] and reflecting the crossover to 
effective single-impurity behaviour (which we emphasise arises whether the system is a 
HF metal or a Kondo insulator). 

2 1 1 1 



3 




Figure 1. T — scaling spectra ttA^D' (ui) and D c (io) vs Cj = lu/ujl for the HCL, 
with r\ = and e c = (dashed), 0.3 (solid). The inset shows the /-spectra on an 
enlarged u> scale; showing that the spectral tails are common, independent of e c . 

At finite temperatures, what one expects for the strong coupling scaling spectra is 
clear: they should now depend universally on a) = u/ojl andT = T/ul- That this arises 
correctly within the present LMA is shown in figure El For a fixed T = 2, the /- and 
c- spectra are shown for progressively increasing interaction strengths U = 5.1,6.1 and 
6.6 with V 2 = 0.2; for e c = 0.3 and rj = (corresponding results for the p-h symmetric 
limit have been obtained in [33]). The inset shows the /-spectra on an absolute scale 
(vs u/U), where the exponential reduction of the a>L-scale with increasing U is clearly 
seen from the change in the width of the resonance. The main figures by contrast 
show the spectra as functions of &, from which the [/-independent scaling collapse is 
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evident; repeating the calculations with different V 2 likewise shows the scaling to be 
independent of V. This behaviour is not of course confined to the chosen T, and figure 
El shows the resultant LMA scaling spectra for a range of T (again for the representative 
e c = 0.3,77 = 0). 




n l 1 l 1 l 1 l 1 l 

-20 -10 10 20 




o 1 1 1 1 1 1 1 1 1 

-20 -10 10 20 

co/co L 

Figure 2. nA a D f (u) and D c (u) vs lo/uj l for the HCL at a fixed f = T/uj l = 2, 
showing scaling collapse with increasing interaction U — 5.1 (dotted), 6.1 (dashed) and 
6.6(solid); for r\ = and e c = 0.3. The inset shows the /-spectra on an absolute scale, 

vs U)/t*. 

Figures El and El show clearly the thermal broadening and ultimate collapse of the 
/-resonance with increasing T; which is naturally accompanied by a redistribution of 
spectral weight leading to infilling of the (u > 0) spectral gap seen in figure ^ for T = 0. 
In fact by T ~ 1 this gap is already obliterated, and the lattice Kondo resonance also 
significantly eroded. This behaviour is typical of the metallic HF state. By contrast, 
corresponding results for the p-h symmetric Kondo insulator (e c = = rj) are shown in 
figure 4 of [33]. In that case the insulating gap at the Fermi level fills up with increasing 
temperature, and the Fermi level D^(u = 0) in particular increases monotonically with 
increasing temperature; in contrast to the the asymmetric HF spectra shown above 
where D*(0) diminishes with T. 
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Figure 3. Thermal evolution of HCL scaling spectra 7rAoD^(o;) and D c (uj) vs 
u> = lo/ujl for e c = 0.3, rj = and temperatures T = 0.2(solid), 0.5 (short-dashed), 
l(double point-dash), 5(dotted), 10(long dash) and 20 (point-dash). Inset shows the 
/-spectra on an enlarged ui scale. 

Two further points regarding figure El should be noted. First, the thermal evolution 
of the /- and c-spectra differ somewhat in terms of the persistence of a pseudogap - the 
/-spectrum shows no sign of the gap by T = 1, while a weak pseudogap structure persists 
in the c-spectrum up to T ~ 5; this reflects the rapid spread of spectral weight caused by 
the meltdown of the sharp resonance in the /-spectra, of which there is no counterpart 
in the c-spectra. Second, the inset to figure El shows the /-spectra on an enlarged 
frequency scale out to u> ~ 100, from which it is seen that the high frequency behaviour 
of the finite-T scaling spectra coincide with that for T = 0. This is physically natural, 
since one expects the dominant influence of temperature to be confined to frequencies 
| a) | < T. The corollary of course is that non-universal frequencies are affected only on 
non- universal, and thus in general physically irrelevant, temperature scales (as shown 
in figure 5 of [33] for the p-h symmetric case). 
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3.1. Scattering rates 

We consider now the scattering rates t -1 (co>;T) that underlie the evolution of the 
conductivity, and are given explicitly in terms of the /-electron self-energy by equation 
1)2. lip . Since the system is a Fermi liquid with £}(cj;T = 0) ~ 0(lu 2 ), at T = 
there is of course no scattering at the Fermi level, r _1 (0;0) = 0. The low frequency 
behaviour of the T = scattering rate can be understood qualitatively by using the 
low-o; expansion of £y (u; 0) (equation ()2.12j) ) and simply neglecting the imaginary part 
Ej(u;;T); leading to 

T- 1 (uo;T = 0)^7i5(u-i* f ) (3.1) 

with e*j = (ef + Sj(0))/l^ 2 the renormalized level. Restoring the small but strictly non- 
vanishing £j(u) ~ iy, T = 0) naturally implies a narrow resonance centred on uj ~ 
instead of a pure 5-function. At finite temperature, we likewise expect the scattering 
rate to increase from zero in the neighbourhood of the Fermi level, reflecting the finite-T 
contribution to S^(u; ~ 0; T); and that this will simultaneously lead to further, thermal 
broadening of the resonance at e^. 

The above picture is corroborated by LMA results as shown in figure displaying 
the u = oj/ujl dependence of r^ 1 (c<j;T) (in units of t* = 1) arising in strong coupling 
for e c = 0.3 and 77 = 0, for a range of temperatures T = T/ul between and 20. In this 
case the renormalized level is found to be ~ 1, precisely where r _1 (co>;T = 0) has a 
narrow resonance. With increasing temperature the resonance is indeed seen to broaden 
and decrease in intensity; and we reiterate that this occurs for temperatures T set by 
the scale ujl - the sole low-energy scale characteristic of the problem in strong coupling. 
Excepting the lowest T we also note that scattering rates in the vicinity of the Fermi 
level are on the order of 0.1 — 1 of the bandwidth t*, values some two or so orders of 
magnitude higher than for conventional clean metals (and indicative of the higher d.c. 
resistivities that are typical of heavy fermion materials [1]). Neither is this behaviour 
confined to a narrow T regime since even for T ^> 1 the scattering rates decay very 
slowly with T; the Fermi level scattering rate for example is readily shown to decay as 
t-\0;T) oc l/ln 2 (T). 

The scattering rates are also related to the /-electron scaling spectra considered 
above. For the Kondo insulating p-h symmetric PAM, it was shown in [33] that the 
dimensionless scattering rate defined as 

—, — = —7 — Tpr = 7/(w T) 3.2 
t(u;T) t(uj;T) 

coincides asymptotically with the /-spectral function, specifically 
1 



~ ttAoD'Cw) (3.3) 
t(uj; 1 ) 

in the regime \Cu\ ^> 1 for any T (the spectral 'tails'), and for all \uj\ for sufficiently large 
T ^> 1. Equation ()3.3p is in fact readily shown to be quite general, and not dependent 
on p-h symmetry. That it holds for HF metals embodied in the asymmetric PAM is 




Figure 4. Thermal evolution of the strong coupling scattering rate T~ 1 (ui;T)/t* vs 
cj/ul (for the HCL) with e c = 0.3, r) = 0; for temperatures T = (solid), 0.5 (dotted), 
1 (short dash), 2 (point-dash), 5 (long dash) and 20 (double point-dash). Inset: the 
same on an enlarged u) scale, showing that the high-(D behaviour coincides with that 
for T = 0. 



illustrated in figure where for e c = 0.3 and r] = the strong coupling f~ 1 (c<j; T) and 
itA D^(uj) vs uj are compared, for T = in the left panel and T = 2 and 10 in the right 
panel. The high-frequency behaviour of the scaling spectrum tiAqD^ (u) is itself known, 
being given (here for rj = explicitly) by [34] 

ttAo^Co;) - ( 1 — 2 + — 2 1 (3.4) 

2^[iln(aM)] 2 + l [iln(aP|)] 2 + 2 5 ; 

with a a pure constant 0(1). These slowly decaying logarithmic tails are evident in 
figure 03 and as mentioned in section 3 embody the connection to effective incoherent 
single-impurity physics on high energy scales. They are independent of the interaction 
U, local hybridization V, underlying conduction band asymmetry e c , and even of the 
lattice type; depending, albeit weakly, only on the /-level asymmetry [34]. 



4. DC transport 

The above discussion of scattering rates leads naturally to consideration of transport; 
beginning with the d.c. limit where (section 2.2) the static conductivity a(0;T) = 
|o"o-F(0;T), with F(uj;T) given for the hypercubic lattice by equation ()2.10j) . In the 
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Figure 5. ttAqD-' (oj) (solid curves) and f _1 (cj;T) (dotted) versus u) (for the HCL 
with e c = 0.3, r\ — 0) at T = in the left panel; and for T = 2 (light curves, solid and 
dotted) and T = 10 (dark curves, solid and dotted) in the right panel. 

strong coupling regime we expect static transport to exhibit universal scaling in terms 
of T = T/ul, and our aim here is to understand its thermal evolution across the full T 
range. Transport on non-universal temperature scales T ~ A (= nV 2 po(— e c )) or ~ £*, 
will be discussed briefly at the end of the section. 

For the p-h symmetric Kondo insulator, LMA results for the T-dependence of the 
scaling resistivity have been considered in [33] (in this case ujl = ZV 2 is equivalently 
the insulating gap scale 'A g '). The T = resistivity is naturally infinite reflecting the 
gapped ground state, the scaling resistivity p(T) = 1/Fhcl(0', T) has an activated form 
p(T) cx exp(a/T) for T <C 1 (with a a pure constant 0{1) and hence a 'transport gap' 
of auij] and p(T) decreases monotonically with increasing T, tending asymptotically 
to incoherent single- impurity scaling behaviour ( [33] and figures 17181 below) . 

For the general case of heavy fermion metals the situation is of course quite different, 
and what one expects in qualitative terms well known [1,2]. The T = resistivity 
vanishes, reflecting the absence of Fermi level scattering and the underlying coherence 
generic to any Bloch decomposable lattice. With increasing temperature p(T) increases 
(initially as ~ T 2 for T < 1 [2,48]), passes through a maximum at ^~niax ^ classic 
signature of HF compounds [1-6] - and decreases thereafter in the strong coupling, 
Kondo lattice regime of interest. Figure |B1 shows LMA results for p(T) vs T for fixed 
e c = 0.3, 77 = 0, and with increasing interaction U = 4.6,5.1,6.1 and 6.6 for V 2 = 0.2. 
The scaling collapse is clearly evident: while the low-energy scale uol itself diminishes 
exponentially on increasing U, universal scaling of p(T) as a function of T — T/ul indeed 
arises in strong coupling, independent of interaction strength (and likewise readily shown 
to be ^-independent on repeating the calculations varying V 2 ). 




Figure 6. The d.c. resistivity p(T) = l/F HCi (0;T) vs. f = T/lu l for e c = 0.3, r/ = 
and four interaction strengths: U — 4.6 (long dash), 5.1 (dotted), 6.1 (short dash) and 
6.6 (solid). Arrows show the inflection points. Inset (a): As in the main figure but on 
a logarithmic T scale, showing the 'log-linear' regime. Inset (b): Scaling resistivity for 
the Bcthe Lattice, for e c = 0.3, T] = 0; the cross shows the T = resistivity of the free 
(V = 0) conduction band. 



This leads us first to comment briefly on the issue of 'the coherence scale', 
characterising the crossover from low-temperature lattice coherent behaviour to 
high-temperature effective single- impurity behaviour. Experimentally, many such 
identifications of the low-energy scale are commonly employed. Some groups use T max 
at which p{T) peaks, others identify the scale via the inflection points (d 2 p(T) / dT 2 = 0, 
marked by arrows in figure EJ), via the leading ~ T 2 behaviour of p(T) at low-T, or via 
the onset of the 'log-linear' regime [26] (shown in inset (a) to figure |B1 and seen in many 
experimental systems [1-6]); the inverse of the T ~ paramagnetic susceptibility, or the 
width of the lattice Kondo resonance, are other possibilities. This leads to what at first 
sight might seem a plethora of low-energy scales. The key point however is that, because 
physical properties in strong coupling scale universally in terms of one low-energy scale, 
all the above definitions of 'the coherence scale' are fundamentally equivalent: all are 
proportional to Ul, and hence to each other — in figure |H1 for example, the inflection 
points in p(T) lie at T = T/ujl = 0.35 and 1.85, and the peak maximum at T = 0.88. 

As for single-particle dynamics and scattering rates considered in section 3, the 
T-dependent scaling resistivity is independent of U or V (as above) but depends 
in general on e c (reflecting the conduction band asymmetry and determining n c via 
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equation (|2.17ojl ) and rj (reflecting the /-level asymmetry). To consider this figured 
shows the resultant scaling resistivities p(T) vs T for rj = 0, and a range of different 
e c = 0,0.1,0.3,0.5 and 0.6, corresponding respectively to conduction band fillings 
n c = 1,0.89,0.68,0.49 and 0.42. The e c = example is the Kondo insulator [33], with 
its characteristic diverging p(T) as T — > 0. The others are all HF metals, and exhibit 
the same qualitative behaviour for all e c — a positive slope for T < T max , going through 
the maximum and then decreasing monotonically for T > T max ; the coherence peak 
itself increasing monotonically with e c , albeit slowly such that T max = T max /ui ~ 0(1) 
for the e c -range shown. Qualitatively similar behaviour is found on varying the f- 
level asymmetry r) for fixed conduction band asymmetry embodied in e c , although 
quantitatively this effect is appreciably less. 




Figure 7. Strong coupling HCL p(T) vs. T for rj = and e c = (solid), 0.1 (dotted), 
0.3 (long dash), 0.5 (double point dash) and 0.6 (short dash) Inset: for e c = 0.5, full 
result (solid) compared to the approximation equation l|4.1|l (dashed) . 

The significant e c -dependence of p(T) = 1/F h cl(0] T) seen in figure[7|for T > T max 
is intuitively natural: the strong coupling Kondo lattice regime corresponds to rif = 1, 
but with variable conduction band filling (n c ) controlled by e c (equation ()2.17ajl ): and 
on decreasing n c (increasing e c ) one expects the static conductivity to diminish and 
hence an increased p(T), as found. To understand the e c -dependence, and in turn 
to enable connection to incoherent effective single-impurity behaviour at high-T, we 
first consider an approximate evaluation of Fhcl(0;T) (equation ()2.10|) ): in which the 
energy dependence of the free conduction band dos g?q(c) = p ( e — e c) is neglected, 
G?o(e) — ^o(0) = Po(— e c) being replaced by its Fermi level value. Employing this 'flat 
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band' approximation in equation (|2.1(J|) (where it enters via the (...) e average) leads to 

F HCZ (0;T)^|[p (-6 c )] 2 ^^^^f(a;;r) = \[p,{-e c )f <f> (4.1) 

expressed as a physically intuitive thermal average of the dimensionless scattering time 
f(u;T) (equation ()3.2)) ). For Kondo insulators this approximation is qualitatively 
inadequate at low-T [33], but as illustrated in figure (inset) it is entirely respectable 
for the HF metals and in particular recovers precisely the high-T asymptotics of p(T). 
As shown in section 3.1, the large Q and/or T dependence of the reduced scattering 
rate f -1 (w;T) coincides with the /-spectral function tiA D^(uj) (equation (|3.3jl ): and 
in section 3 (see also [34]) the latter were shown to have common spectral tails, 
independently of e c . This suggests that the primary effect of e c seen in figure [7| for 
p(T) = 1/Fhcl(0;T) is contained in the [po(—£ c )} 2 of equation (|4.1jl . 

That this is so is seen in figure |H1 where the results of figure [7| are now shown as 
p'(T) vs T, where 

'(T) - ^ Po( ~ ec)]2 U 2) 

For T > 5 or so in practice, p'(T) is seen in particular to be independent of the 
conduction band filling embodied in e c ; including we note the Kondo insulator, whose 
'high' temperature resistivity is thus seen to be that of a regular heavy fermion metal. 
Indeed as readily demonstrated, and evident in part from the above discussion, the 
behaviour seen in figure |H1 is barely dependent on the details (e-dependence) of the host 
bandstructure embodied in po( e )- 

The obvious final question here concerns the high-T form of p'(T) for the PAM. To 
that end we consider the Anderson single-impurity model (AIM), with p imp (T) denoting 
as usual the change of resistivity due to addition of the impurity to the non-interacting 
host, and p[ mp (T) = Pi mp (T) / p imp (0) . This is given by [2] (cf equations ()4.1|) . (j4.2|) ) 

/ ( r p\ = / duJ —a r im P 4.3 

with the impurity scattering rate f^L(u;T) = irA D imp (Lu;T); where D imp (w;T) is 
the impurity spectral function such that 7rA _D imp (0; 0) = 1 follows from the Friedel 
sum rule [2] in the singly occupied, strong coupling Kondo regime of the AIM. The 
LMA scaling resistivity p' imp (T) vs T is also shown in figure |Hl where T = T/ujk and 
ujk = Z imp V 2 is the AIM Kondo scale (with Z imp the impurity quasiparticle weight). 
From this it is seen that the high-T scaling behaviour of p'(T) for the PAM is precisely 
that of the AIM; in particular the leading T>1 behaviour of the LMA p'(T) is readily 
shown analytically to be given by p'(T) ~ 3tt 2 / (16 ln 2 (T)), which is exact in the Kondo 
limit of the impurity model [2] . This reflects again the crossover in the strong coupling 
PAM from low-temperature lattice coherent behaviour to incoherent effective single- 
impurity scaling physics, here in the context of d.c. transport. As for its counterpart 
in the case of single-particle dynamics [34], we point out (a) that since this connection 




Figure 8. Scaling resistivities p'(T) vs. T for rj = and e c — (solid), 0.3 (long dash), 
0.5 (dotted) and 0.6 (point dash); c.f. figure [7| p[ mp (T) for the single-impurity AIM 
(equation 14.3J0 is also shown (double point dash). For inset, see text. 



is established from scaling considerations it is entirely independent of how the scales 
ojl and u K for the two distinct models (PAM and AIM) depend on the underlying 
bare/material parameters of the respective problems; and (b) the fact that it arises in 
the T/ljl scaling regime precludes a two-scale description of the crossover from lattice- 
coherent to incoherent effective single- impurity physics. 

Our focus above has naturally been on the strong coupling, Kondo lattice regime. 
We now look briefly at d.c. transport on non-universal scales. What one expects here 
is that when the temperature is a not insignificant fraction of the hybridization Ao or 
bandwidth scale £*(= 1), Kondo screening will be washed out, and hence p(T) should 
cross over from the logarithmically decreasing single-impurity form at T ^> 1 (figures [7| 
and IE} to conventional metallic behaviour dp{T)/dT > at non-universal temperatures; 
and thus as such must go through a minimum. That this indeed happens can be seen 
in the inset to figure |E1 where we show p{T) = 1/Fhcl(0;T) vs T = T/ujl for e c = 0.3, 
i] = 0, V 2 = 0.2 and U = 4.1 (solid line) and 5.6 (point-dash). For the lower U example, 
a minimum is seen at T = T/u L ~ 300, which corresponds in 'absolute' units (£*) to 
a temperature T ~ 0.2 — an appreciable fraction of the hybridization Ao — 0.3. The 
corresponding minimum does of course exist for the higher U, but is pushed beyond 
T = 10 3 (and to concomitantly lower vales of p(T)); and p(T) in this case lies on the 
universal scaling curve throughout the T-range shown in figure |H1 

A final point is worth noting here. For T < 10 2 , the p(T) vs T = T/ul for the 
two Vs shown in figure |E1 (inset) are in essence coincident; each lies on the universal 
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scaling curve. What distinguishes different interaction strengths is of course the location 
of the minimum, occurring as it does on non-universal temperature scales. No real 
HF material is however in the universal scaling regime 'for ever' — with increasing T 
the scaling regime will be exited sooner or later. And the temperature for which the 
experimental p(T) is a minimum (once phonon contributions have been subtracted out) 
can provide valuable information on the interaction strength, as we shall see in action 
in the following paper. 

4-1. Bethe lattice 

We have considered almost exclusively the hypercubic lattice, for the obvious reason that 
its one-particle Bloch states ultimately underlie the low-temperature lattice coherence 
of the interacting problem. For the Bethe lattice, the strong coupling scaling resistivity 
p(T) = 1/F BL (0;T) (with F BL (u; T) from equation (|2.18|l ) is shown vs T in inset (b) to 
figure^ for e c = 0.3, n = 0. In contrast to its counterpart for the HCL shown in the main 
figure, the T = resistivity is non- vanishing (given by equation (|2.19j) ). reflecting the 
absence of coherent Bloch states for the BL. Further, the high-T asymptote of the BL 
p(T) in the scaling regime is likewise non- zero; being given by the T = value of the free 
(V = 0) conduction band resistivity, namely l/[po(— e c )] 2 = 7r 2 /[4(l — e 2 )] as marked by 
a cross in figure El inset (and arising for the same physical reasons discussed for Kondo 
insulators in [33]). The qualitative contrast between p(T) for the canonically Bloch 
decomposable HCL, and that for the BL, illustrates why the latter — more specifically 
the associated 'joint density of states' type formula equation (j2.18|) for F(0; T) that is 
not uncommonly employed in the literature — gives a poor caricature of d.c. transport 
for HF metals in which the lattice coherence is of central importance. 

5. Optical conductivity 

We turn now to the optical conductivity a(u; T) = ^a F(u}; T) (with Fhcl(&] T) given 
by equation (|2.1()j) ). In the strong coupling Kondo lattice regime Fhcl(u; T) is of course 
independent of U and V 2 , and a universal function of u — u/u L and T = T/ul for fixed 
e c and rj. 

LMA results for Fhcl{u;T) are shown in figure EH for e c = 0.3 and 77 = 0. The 
right panel shows the thermal evolution of the optical conductivity (on a linear d)-scale) 
for temperatures T = 0, 0.5, 1, 2, 5 and 10; while the left panel (on a log-log scale) shows 
the behaviour for a lower range of temperatures up to T = 0.5. The latter in particular 
illustrates the thermal evolution of the optical Drude peak, which at T = consists of an 
uj = 5-function given by equations (2.14) (with net weight oc uj^ in strong coupling). On 
increasing T from the Drude peak naturally broadens, and is well fit by a Lorentzian 
up to its half-width or so, after which it decays more slowly in Co. At the lowest T 
shown the Drude peak is well separated from the 'optical edge' in Fhcl^iT) seen 
at Co ~ 2 (although we add that Fhcl{u;T) is strictly non-zero for all uj), and with 
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Figure 9. Fhcl(w,T) vs. ui = lu/lol in the Kondo lattice scaling regime for a range 
of temperatures T — T/u>l', for e c = 0.3 and 77 = 0. Left panel: On a log-log scale, for 
f = 0.02 (solid), 0.1 (dotted), 0.2 (short dash), 0.3 (point dash), 0.4 (long dash) and 
0.5 (double point dash), showing the evolution of the Drude peak. Right panel: On 
a linear scale, for T = (solid), 0.5 (dotted), 1 (short dash), 2 (point dash), 5 (long 
dash) and 10 (double point dash). 



increasing T is seen to persist as an essentially separate entity up to T ~ 0.1 or so; after 
which it is progressively destroyed as expected, merging into an optical pseudogap in 
the neighbourhood of £0 ~ 1 — 2, which is reasonably well filled up by T ~ 0.5 and all 
but gone by T ~ 2 (see figure El right panel). Similar behaviour is naturally found on 
varying e c and/or 77. Figure [TU1 shows in particular the influence of e c (varying conduction 
band filling) on the optical pseudogap for a fixed temperature T = 0.2, from which it is 
seen that the pseudogap becomes shallower with increasing e c . 

The above behaviour should be compared to the p-h symmetric Kondo insulator 
(KI) e c = = 77 considered in [33]. In that case the T = optical conductivity is 
characterized by an indirect gap A ind = 2ZV 2 — 2u>l, and there is of course no T = 
Drude peak. Instead a Drude-like peak in the optical conductivity actually builds up on 
initially increasing T from zero (see figure 15 of [33]), before being thermally broadened 
and subsumed into the optical pseudogap. For T < 1 or so the low-frequency optics of 
the KI are thus very different from those of the HF metal, as expected. But for T > 1 
the optical behaviour of the two is qualitatively similar as shown by comparison of figure 
E3 (right panel) and its counterpart for the KI, figure 14 of [33]. This too is physically 
natural, since the infilling of the indirect optical gap on temperature scales T ~ 1 means 
that the KI behaves to all intents and purposes as a HF metal; as seen also in figures U\ 
or|Hlfor the static transport. 

A second point should be emphasised here, obvious though it is from the preceding 
discussion: whether for HF metals or Kondo insulators, it is the low-energy scale 
u>l = ZV 2 that sets the intrinsic scale for both the w-dependence of the low-energy 
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Figure 10. Fhcl(u;T) vs. u) for fixed T — 0.2 in the Kondo lattice scaling regime 
for r\ — and e c =0.1 (solid), 0.3 (dotted), 0.5 (short dash), 0.6 (dot dash) and 0.8 
(long dash). Inset: shown on a linear scale. 



optical conductivity and its thermal evolution. And in strong coupling that scale is 
wholly distinct from the optical direct gap, A dir . The latter arises at its simplest 
in the commonly employed renormalized band picture (see e.g. [2]), as the minimum 
direct gap for which optical transitions are allowed. In this effective single-particle 
description the imaginary part of the /-electron self-energy - and hence all scattering 
- is neglected entirely, and the corresponding real part £j (a>; 0) is replaced by its 
leading low-a; behaviour equation (j2.12|) (as also inherent to a slave boson mean-field 
approximation [2]). The two branches of the renormalized bandstructure, denoted by 
ct>±(e) with e = ek, then follow from the zeros of [G c {e^] uj)}^ 1 = [7(0;) — (see equation 
(|2.7Jl ) with the approximate (7(0;) =) Jr(uj) ~ uj — e c — ZV 2 [uj — Ze*^ 1 from equation 



(|2.56j) : and the resultant e-dependent direct gap A^(e) = [cJ+(e) — ^-(e)] is given by 
A d (e) = [(e + e c -Ze* f ) 2 + AZV 2 ] 1/2 with e* f = e f + Ef(0;0) the usual renormalized level. 
The minimum direct gap, A^r, occurs for e + e c = Ze*j (~ in strong coupling) and is 
thus 



Mir 



2\fZV. (5.1) 



The corresponding result for the optical conductivity Fhcl(^',T) is readily determined 
from equation ()2.10|) . Denoted by F (uj; T) it is given for T = (and all uj > 0) by 



F (cu; 0) 



9(uj - A dir ) A^ ir 



d o ( 4 + V^ 2 - A dir )+d c [e}-Juj 2 - A 2 dir 



(5.2) 



with d^iuj) = p (uj — e c ) the free conduction band dos and 9{x) the unit step function; 
and is thus non-zero only for frequencies uj > A^ ir above the direct gap (which result is 
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also readily shown to hold for all temperatures). 

Two points should be noted here. First that the low-energy scale ujl = ZV 2 
intrinsic to HFs or KIs is qualitatively distinct from the direct gap Adi r . In fact since 
Adir/^i oc l/y/uJ£ it follows that in strong coupling where the quasiparticle weight Z 
and hence ujl becomes exponentially small, optics on the direct gap scale do not even lie 
in the to = uj/ujl scaling regime; although neither do they occur on truly non- universal 
scales (because A dir oc \fz) and in that sense belong to the 'low-frequency' optical 
spectrum. Second, we emphasise the inherent naivete of interpreting optics in terms 
of renormalized single-particle interband transitions: it is scattering due to electron 
interactions that generates all the optical density below the direct gap scale. Failure to 
include such, as in a renormalized band picture — and regardless of how sophisticated 
the underlying band structure employed in practice — inevitably leads to a qualitatively 
inadequate description of optics (as illustrated explicitly in figure [TT] below). Neither 
is this situation ameliorated in materials application by the introduction of ad hoc uo- 
dependent broadening factors, for that simply avoids the basic underlying physics. 



Figure 11. Fhcl(u; T) vs Cj = uo/lol on a log scale for e c = 0.3, rj = with U = 6.6 
and V 2 = 0.2; and temperatures T — T/lol — (solid), 0.5 (dotted), 1 (short dash), 
2 (point dash), 5 (long dash) and 10 (double point dash). Comparison is also made 
to the renormalized band picture equation (|5.2|l (thick solid line). Inset: renormalized 
bandstructure w±(e)/wi versus the free (V = 0) conduction band energies e = e^. 



LMA results for optics on all frequency scales are given in figure for e c = 
0.3, rj = 0; where F HC l{u; T) is shown vs u — lu/ul on a log scale spanning five orders 
of magnitude, for the same range of temperatures T = T/ujl employed in figure El (right 
panel). To encompass all u including non-universal energies, the bare parameters U and 
V 2 must of course be specified, U = 6.6 and V 2 = 0.2 here being chosen for illustration; 
although note that the optical conductivity as a function of to = uj/ujl remains 'universal' 




3 



10 



4 
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Figure 12. Fhcl(w,T) vs. ui = uo/lol for the same parameters as in figure ITT1 and 
temperatures T = T/ljl from to T — 500 ~ §Adi r as marked on the figure. Inset: 
same on a linear scale. 

(independent of U or V 2 ) up to large but finite values of £u determined by the particular 
U and V 2 chosen, in this example u ~ 300 — 400 (a directly analogous situation for 
the resistivity p(T) is shown in the inset to figure The inset to figure HP shows the 
renormalized bandstructure o)±(e) = uj±{e)/ujL vs the free conduction band energy e; 
determined as above from solution of Jr(uj) = e (with 7r(u;) the full Re(7(u;))). This 
enables the notional direct gap to be determined, Adir — 750ul here — well separated 
from the low-energy coherence scale ujl in strong coupling — and indeed seen to occur 
for e ~ — e c = —0.3. 

The essential points from figure ^2 are clear. As expected and well known (see 
e.g. [11, 12]), significant optical absorption occurs in the vicinity of the direct gap; 
strongly broadened to low-energies due to electron interactions as above, and all but 
'dead' on non-universal energy scales (e.g. the hybridization A = nV 2 p (— e c ) ~ 
7 x 10 3 ujl for the chosen bare parameters). Regarding the thermal evolution of the 
optical conductivity note also that temperatures on the order of a few multiples of the 
coherence scale uo^ — which control the thermal evolution of the low-energy optics - 
have essentially no effect on frequencies of the order of the direct gap, reflecting the 
clean separation between ul and Adi r characteristic of strong coupling. As a corollary 
the direct gap should be thermally eroded only for T ~ O(Adir)', as indeed seen in figure 
H~2*l where (for the same parameters as figure ITT|l the thermal evolution of Fhcl{^',T) 
is shown for temperatures up to T = 500 ~ |Adi r . Significant thermal erosion sets in 
by about T/A<ar ~ 0.2 or so, and is well developed by the highest temperature shown. 
The clear scale separation between ojl and Adi r will not however be captured properly if 
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one is restricted to relatively low interactions and high temperatures as e.g. in quantum 
Monte Carlo [16,18], or from theories in which the quasiparticle weight Z is algebraically 
rather than exponentially small in the interaction strength, such as iterated perturbation 
theory [22,23]. 

6. Conclusion 

We have considered here the periodic Anderson lattice, the canonical model for 
understanding heavy fermion metals, Kondo insulators, intermediate valence and related 
materials. Optical conductivities, d.c. transport and single-particle dynamics of the 
paramagnetic phase have been investigated, using the local moment approach within a 
DMFT framework. For obvious physical reasons our main focus has been the strongly 
correlated Kondo lattice regime, where we find the problem to be characterised by a 
single, exponentially small coherence scale ul] in terms of which the frequency and 
temperature dependence of physical properties scale — being universally dependent on 
uj = uj/ujl and/or T = T/ujl regardless of the interaction or hybridization strengths. 
All relevant energy/temperature scales are handled by the theory, from the low-energy 
coherent Fermi liquid domain out to large (and in the strict scaling limit arbitrarily large) 
multiples of lvl where incoherent many-body scattering dominates the physics; followed 
by the crossover out of the scaling regime to non-universal, high energy/temperature 
scales dictated by 'bare' model/material parameters. And while our emphasis has 
been on strong correlations we add that all interaction strengths from weak to strong 
coupling are encompassed by the LMA [34], such that intermediate valence behaviour 
in particular can also be addressed. 

The first question posed in the Introduction nonetheless remains: to what extent 
does the model, and our theory for it, capture experiment? We turn to that in the 
following paper where direct comparison of theory and experiment is made for three 
heavy fermion materials and a classic intermediate valence compound. 
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